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We study the impact of impurities on turbulence driven intrinsic rotation (via residual 
stress) in the context of the quasi-linear theory. A two-fluid formulation for main and 
impurity ions is employed to study ion temperature gradient modes in sheared slab 
geometry modified by the presence of impurities. An effective form of the parallel 
Reynolds stress is derived in the center of mass frame of a coupled main ion-impurity 
system. Analyses show that the contents and the radial prohle of impurities have a 
strong influence on the residual stress. In particular, an impurity profile aligned with 
that of main ions is shown to cause a considerable reduction of the residual stress, 
which may lead to the reduction of turbulence driven intrinsic rotation. 
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I. INTRODUCTION 


The presence of a certain amount of impurities is, in some sense, an inevitable conse¬ 
quence of tokamak plasma operation. The contents and the radial prohle of the dominant 
impurity species depend on the details of discharge conditions, plasma-wall interaction, and 
the physical process of impurity transport. In tokamaks, the material constituting plasma 
facing components (PFCs), for instance, carbon in DIII-EP and KSTARP, or tungsten in 
ASDEX-lP, JETHI and ITERP, is usually observed as a dominant impurity species. In addi¬ 
tion to impurities originating from the surrounding PFCs, a considerable amount of helium 
ash will be present in ITER or fusion reactors as a result of fusion reactiorP^. 

Studies of impurity transport in tokamaks have usually focused on the determination of 
impurity prohles (e.g. the impurity peaking factor) from either neoclassical®^ or turbulent 
mechanism^n^. The main goal of this endeavour is to predict the degree of impurity 
accumulation for a given set of discharge parameters to avoid deleterious effect from impurity 
accumulation. Turbulent transport of impurities, which has been a subject of an extensive 
study in fusion plasma physics for decade^^^®^ is rekindling interest in this regard since 
recent experimental evidence has shown the importance of the turbulent transport process 
in determining the core impurity prohld^^. 

When impurity concentration is sufficiently low, one may treat the impurities as passive 
reactants to background turbulence. In this case, the so-called trace approximation has been 
used to study turbulent impurity transport !^^*^^l If impurity concentration is considerable, 
however, they start to affect microinstabilitj^^®^ and the trace approximation becomes 
invalid, as pointed out by Fiilop and NordmarP®. Most notably, the alignment of a impurity 
profile with that of main ions has been known to yield a significant change of stability of 
the ion temperature gradient (ITG) mode, sometimes leading to the destabilization of an 
independent impurity drift mod#^*^ 

Our interest in this paper lies in the study of the impact of impurities on turbulent mo¬ 
mentum transport, in particular, the generation of intrinsic rotation via residual stresd^^®. 
We note that few results are available in this aspect, given a lot of published papers on 
turbulent impurity transport in tokamaks. Intrinsic rotation in magnetic fusion plasmas 
is thought to be generated through the conversion of radial inhomogeneity into (= 
spectrally averaged wavenumber in parallel direction) asymmetry via residual stress. This 
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conversion process requires some symmetry breaking mechanismP®!! early paper cal¬ 
culated the particle flux in the presence of combined ITG, impurity and parallel velocity 
gradients, but not addressed the turbulent residual stress giving rise to intrinsic rotatiorP^. 

The purpose of this study is to provide a simple physics insight into which how turbulence 
driven residual stress is related to the characteristics of impurities. Intuitively, one may ex¬ 
pect that the change of stability and characteristics of an unstable mode structure in the 
presence of a considerable amount of impurities will alter the turbulent residual stress. This 
effect may change the amount of intrinsic torque anticipated in toakamak experiments. Since 
intrinsic rotation is envisioned to provide necessary plasma rotation in reactor-relevant toka- 
mak experiments, including ITER, because of the limited capability of driving the necessary 
torque in a reactor by neutral beam injection, it is of importance to conduct a sufficiently 
detailed study of this effect. 

To address the problem mentioned above, we employ a set two fluid equations for main 
and impurity ions. By treating impurities as a separate species satisfying impurity fluid 
equations, we can incorporate the effects of impurities on the characteristics of the ITG 
mode self-consistently. For simplicity, we adopt the sheared slab geometry. By avoiding the 
complication in analysis due to the sophisticated toroidal geometry, we focus on the basic 
physics of the impurity effect on residual stress in the presence of impurity-modified ITG 
turbulence. Analyses show that the relative importance of residual stress to the diffusive 
momentum flux is reduced when the gradient of impurity and main ions has the same sign. 
This indicates a possible decrease of intrinsic rotation when a considerable amount of core- 
peaked impurity is present. 

The remainder of the paper is organized as follows. In Sec. II, we describe the basic 
formulation. A derivation is given of a linear dispersion relation and an eigenmode equation 
for ITG modes using two fluid equations for ions and impurities. In Sec. Ill, we perform a 
numerical analysis to evaluate eigenmode structures of unstable modes. A detailed study 
of the change of eigenmode characteristics, such as the variation of the mode shift off the 
rational surface, is made in this section. Section IV is devoted to the quasi-linear calculations 
of momentum flux. After introducing some underlying assumptions made in this study, we 
calculate turbulent parallel Reynolds stress (both diffusive and residual stresses) in the 
context of the quasi-linear analysis. A particular emphasis is placed on the evaluation of the 
ratio of residual stress to diffusive one, which ameliorates the limitation of the quasi-linear 
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theory. Possible implications of the results will be discussed in this section. We conclude 
this paper in Sec. V with a brief summary of main results and some discussions. 


II. FORMULATION 


We begin with a set of two fluid equations for main and impurity ions consisting of the 
conservation of density (n^), parallel momentum density {mjnjV\\j), and pressure Pj, 


dnj 

dt 


+ V ■ {nj^rj) = 0, 


TTljU^- 


dv 


id 


dP^ 

dt 


dt 


+ V||Pj + njejV||0 = 0, 


+ rP,'V-V|| =0, 


( 1 ) 


where the subscript j = i (Z) denotes the main (impurity) ions, 0 is the electrostatic 
potential fluctuation, and P is the adiabatic index. The fluid velocity is decomposed into 
the parallel and the perpendicular components, vj = vyj + v±j with v±j = + Vpj, 

where v^j, v^j, and Vpj are the E x B, the diamagnetic, and the polarization drift of the 
species j, respectively. We employ the sheared slab geometry with a background magnetic 
field in the vicinity of a reference surface (ro), Bq = Bq{z + x/Lgy). Here, Lg = qR/rq' is 
the magnetic shear scale length with q the safety factor, R (r) the major (minor) radius, 
and the prime denotes a derivative with respect to x = r — vq. 

Linearization of Eq. (1) can be done straightforwardly giving rise to 


df {jii — V^0) + Lei(l + iPiV^) Vy0 + V ||U||i — 0, (2a) 

dfv\\i - Vo\\iVy(t) + V||0 + TiV\\Pi = 0, (2b) 

dfpi + Lei{l + r]i)Vy4> + rV||U||i = 0, (2c) 

dfinz - |vi0) + L,z{l + -^KzVDVyC^ + = 0, (2d) 

Z 7 ~ 

dfv\\z - Vo\\z^y(t) + 

dfpz + L^zi^ + r]z)'^y(p + rV||U||z = 0, (2f) 


where df = {dt + xV^f^y) with U0g the normalized equilibrium E x B flow shear. Various 
quantities in Eq. (2) are normalized as a; = x/p^, y = y/ps, z = zjLnei t = t/{Ln^/cg), 0 = 
{e(i)/Tf){LnJps), rii = {nii/nio){LnJPs), nz = {nzi/nzo){LnJ ps), U|| = (ui||i/cJ(L„^/pJ, 
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Vz\\ = ivz\\l/Cs)iLnJps), = {LnJCs){dVEo/dx), = {LnJCs){dVi 0 \\/dr), = 

{Lr,Jcs){dVzo\\/dr), pi = {Pii/Pio){LnJPs), Pz = {Pzi/Pzo){LnJps) with L-^ = -Vno/no 
the equilibrium density scale length, Cg = -s/Teo/rrii (Teo: equilibrium electron temperature), 
and ps = Cg/fli the ion-acoustic Larmor radius. Also, various dimensionless parameters are 
defined as /i = pi = L^/Lxi with L^l = — VToj/To* the equilibrium ion temperature 

scale length, pz = L^z/Ltz, Ki = Ti{l + pi) with n = Tio/Teo, Kz = rz(l pz) with 
Tz Tzo/TeOi Lei aud Lg^z Lyi^jL^^. 

We consider perturbations of the form / = fk{x)e'x.p{ikyy + ik^z — iojt), where u is 
normalized to Cg/Ln- Then, one can calculate ion and impurity density fluctuations from 
Eqs. (2a) to (2f), 


'^ik — 


CJ 


{CJ “1“ kyL^-i^i) ^J_ kyL^ 


k\\Cj 


— kuTiV 


, kykwLfi-iKi 

h - kyVZ + 


u 




and 


nzk = 


bJ 




kyLgzKz \ , 2 


ki — k„L 


yl-eZ .2 


k\\(l) 


- mrz/p)T 


k\\ — kyVLz + 


kyk^^ Lgz z 


bJ 


where u = lo — kyxV^Q is the normalized frequency of the fluctuation accounting for the 
Doppler shift. 

Ions, impurities, and electrons are coupled through the quasi-neutrality condition, Ugk = 
(1 — fz)nik + fznzk where fz = Znzolngo is the fraction of impurity under consideration. 
Assuming the adiabatic electron response, Ugk = (j>k, one can derive, after a little algebra, 
an eigenvalue equation in (pk, 


d‘^(t)k 

dx"^ 


+ U{k\\,ky,bj)(l)k 


0 , 


( 3 ) 


where 

U{k\\,ky,u) 


= -k, + 


k,, — U! 


+ 


(1 — fz)d! 


^ au! + kyA {au + kyA){u‘^ — TikuT) 


k-,.k?,LpjKj 

ki-kyhV'..+ " II/' 


bJ 


fzdj 


{aul + kyA) bj"^ — {tz/ p)kuT 


— /cy — kyk^^V^^^z P 


kykjLezKz 

pCj 


( 4 ) 


with the definitions a = 1 - fz + {p/Z)fz and A = (1 - fz)LeiKi + fz{p/Z‘^)LezKz- 
Equation (|^ represents a complete potential function for the ITG eigenmode modified by 
the presence of impurities whose fraction is fz- Following Singh et. we assume that 


4’ki 
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the mode frequency is much higher than the equilibrium shearing rate and the ion acoustic 
wave frequency, i.e., 3> /cyrjF and u 3> kyV^Q. We further assume that the equilibrium 

ion and impurity velocity prohles are equal, Vo\\i = Vo\\z = l^o||- Under these assumptions, 
Eq. Q is reduced to 


U{ky, u, x) = Ao + Aix + A 2 x‘^, 


(5) 


where 



f^uj + kyB k‘^s^ 

aoj + kyA xA ’ 


with the dehnitions /? = 1 — /z + [Zj[i)fz and B = {1 — fz)LeiKi + {fz/fi)LezKz- Equa¬ 
tion (5) with the coefficients Aq,Ai, A 2 is a simplihed form of the impurity-modihed potential 
function for the ITG mode with the frequency uj and the mode number ky. The impurity 
effect on mode characteristics are now contained in the parameters a, A and B in this 
formulation. 

Solutions to Eq. ([^ are well-known Hermite polynomials. To make an analytical progress, 
we consider the most dominant zeroth order Hermite polynomial. Then, the eigenfunction 
becomes 



where 



(7) 


is the shift of the mode from the mode rational surface. 




( 8 ) 


is the width of the mode, and 



2 " 


( 9 ) 
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FIG. 1: A schematic diagram illustrating the wall-peaked {L^z < 0, red dotted) and the 
core-peaked {L^z > 0, green dashed) impurity profiles. Solid blue line represents an 

equilibrium plasma density prohle. 


We plug Eq. ([^ into Eq. ([^ and evaluate the resulting relation, 

A2 




which converts into the following dispersion relation, 




_ ^2 , K ^ _ -ji _I _ V' 

^ au + kyA 4D{au + kyAy \kyS "" 


U 


( 10 ) 


where D = {j5oj + kyB)/{auj + kyA). One can rearrange Eq. (10) to obtain 


cj (1 -|- cxky) -|- ky(jj{^ —1 -|- i'\/Ds(y.) -|- i'\/DkySA 


4D{aoj -|- kyA) 


i— 

\kyS 




It is an easy task to show that Eq. ( |TT| ) recovers the dispersion relation without impuritie^ 
when fz = 0,C( = l,A = Ki and D = 1. 

Equation ( pAj ) is the desired equation to be analysed in this paper. It has been known 
that the radial profile of an impurity has a significant influence on ITG stabilitj^^®^^. In 
the case of a wall-peaked impurity prohle schematically shown in Fig. 1 (red dotted line), the 
ITG mode is destabilized while the core-peaked one (green dashed line in Fig. 1), stabilizes 
it. Therefore, before performing a detailed numerical analysis, it is instructive to examine 
if Eq. (0 recovers these known results. To make the analysis simple, we hrst consider 
Eq. ( pT| in the slow mode (he., the long wavelength) limilP^I^ which is characterized by the 
relation \kyA\ 1. We further assume that = 0, and s <C 1. Then, one obtains 
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a purely growing mode with the growth rate, 




ikysA= ikyS 


= ikyS 


(1 — fz)LeiKi + fZ^LezKz 

Ki - (/f, - IjA'z 


( 12 ) 


where the relation L^i = (1 — fzLez) /(I ~ fz) is used from the quasi-neutrality condition. 


From Eq. (12), one can see that the mode becomes stabilized when L^z > 0 (he., the core- 
peaked case), while it is destabilized when L^z < 0 (he., the wall-peaked case). In slab 
ITG modes, this can be interpreted by an effective reduction or enhancement of negative 
compressibility due to impurity inertia. We now consider the other limit where the relation 
|1 — k‘^A\ < s 1 holds, he., the fast mode limit. In this limit, the fastest growing mode 
occurs at ky = A~^, giving rise to the dispersion relation. 


Again, the mode becomes stabilized (destabilized) when Lez > 0 [L^z < 0) due to the 
decrease (increase) of A = (1 — fz)Lf,iKi + {fi/Z‘^)fzL(,zKz. 


III. EIGENMODE ANALYSIS 

Figure [^a) and 2(b) show the real frequency (ur) and the growth rate ( 7 ), respectively, 
as a function of Lez = Ln^/Ln^ for two representative impurity ions in tokamaks, He^’*' in 
blue and C®"*" in red. To produce Fig. [^we solved Eq. ( |TT| ) numerically using the parameters, 
Ti = Tz = 1, fz = 0.3, ky = 0.5, Tji = 2, = —0.01, and = —0.1. The wavenumber 

ky = 0.5 corresponds to the most unstable mode giving rise to the maximum growth rate in 
this case. The equilibrium E x B shear = —0.01 roughly corresponds to Er ~ —5 kV if 
we assume Bt = 1 Tesla and Te ~ 3 keV. This is approximately ~ 1/4 of the radial electric 
held well in typical H-mode discharge^^^^. Deuterium is used as the main ion species in these 
calculations and throughout the paper. Figure 2 basically conhrms the analytic estimation 
of the impact of Lf.z on ITG stability presented in Sec. II: stabilization (destabilization) of 
the ITG mode by a wall-peaked (core-peaked) impurity prohle. The lighter impurity species 
(He^^) shows a little higher real frequency and growth rate compared to a heavier one, but 
the deviation is insignihcant as long as the fz value is hxed. 

The impact of impurity fraction, fz-, on the eigenfrequncy is shown in Figs, [^a) and 3(b). 
In these hgures, we consider two different impurity prohles: L^z = 1 (solid) and L^z = — 1 













FIG. 2: The real frequency {cOr) (a) and the growth rate ( 7 ) of the most unstable 
eigenmode corresponding to ky = 0.5 as a function of the ratio between the equilibrium 
scale length of the electron (L„e = —"n-o/Vno) to the impurity {Lnz = —'^zo/Vn^o) density. 

He^’*' (blue) and (red) are considered as representative impurities in tokamaks. 

Deuterium is assumed to be the main ion species. Plasma parameters being used are 
Ti = Tz = 1, fz = 0.3, rji = 2, = -0.01, and GJn = -0.1. 

(dotted). Other parameters to produce Fig. |^are same as those of Fig. The stabilization 
(Lez = 1) or destabilization (Lez = —1) of the ITG mode by impurities becomes stronger 
as fz increases. An interesting observation is that there is asymmetry in the stabilization 
and destabilization by impurities for a hxed fz value, as can be seen in Fig. [^b). Thus, a 
heavy impurity with a core-peaked prohle has a stronger stabilizing influence compared to 
the lighter impurity, while the destabilization effect due to it with a wall-peaked prohle is 
weaker than that of the latter. 

Of particular interest in this paper is the strength of the mode shift (Xq) given in Eq. ([^. 
This is because Xq driven by the E x B shear is an essential quantity in the calculation of 
residual stres^^ ^^*^’^ * ^^* , hence the intrinsic rotation, in the quasi-linear theory. Figure [^a) 
shows |Xo| vs. Ln^/Lnz for G®+ when = 0.1 (solid) and 0.3 (dotted). In general, |Xo| for 
a given value of decreases as the local impurity prohle changes from wall- to core-peaked 
one, regardless of fz- Thus, one may expect that the ITG-driven residual stress will increase 
as the impurity prohle changes from the wall to core-peaked one. The trend of |Xo| reduction 
shows a remarkable diherence depending on the impurity contents. When fz = 0.1, |Xo| 
shows an almost linear behaviour in the whole range of conceivable Ln^/L^z values, as can 
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FIG. 3: The real frequency (ur) (a) and the growth rate ( 7 ) of the most unstable 
eigenmode as a function of impurity contents (/z). Solid and dotted lines represent when 
Lez = 1 and — 1 , respectively. He^"*" (blue) and C®’*' (green) are considered as representative 
impurities. The same plasma parameters (except for fz) as those in Fig. 2 have been used. 

be seen in Fig. 4(a). When fz = 0.3, however, |Xo| does not change much when negative 
Lrie/Lnz < 0 , while it starts to decrease very rapidly as soon as Ln^/Ln^ becomes positive. 
For example, the reduction in |Xo| as Ln^/Ln^ changes from 0 to 1 is ~ 25% and ~ 9% 
when fz = 0.3 and 0.1, respectively. 

Figure l^b) shows the effect of fz on |Xo| when = 1 (solid) and —1 (dotted), 

respectively. Unlike the growth rate, |Xo| always decreases as fz increases, regardless of 
the impurity prohle shape. The decrease of Xq can be seen more clearly in the core-peaked 
impurity prohle which exhibits a signihcant reduction of |Xo| when fz increases. The impu¬ 
rity effect on Xq is not considerable for a wall-peaked prohle with a moderate dilution of a 
plasma, he.,, when fz < 0 . 2 . A core-peaked prohle, however, will result in large reduction 
of Xq, which eventually leads to a signihcant decrease of residual stress. This is a subject 
for the next section. 


IV. MOMENTUM FLUX AND RESIDUAL STRESS 

In this section, we calculate the parallel Reynolds stress in the presence of a small fraction 
of impurity ions. To begin, we add parallel force balance equations for ions and impurities, 
Eqs. (2b) and (2e), respectively, and take an ensemble average of the resulting equation to 
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FIG. 4: The shift (Xq) of the most unstable ITG mode {ky = 0.5) off the rational surface 
as a function of Lez = Ln^/Ln^ (a) and fz (b). Blue solid and red dotted lines in (a) 
represent when fz = 0.1 and 0.3, respectively. In (b), blue solid and red dotted lines show 
when Lf>z = 1 and —1, respectively. G®"*" impurity is used with plasma parameters 
Ti = Tz = 1, r]i = 2, VfjQ = -0.01, and GJ|| = -0.1. 


obtain 


^(^l> 


dt 


dr 


„(i) mznzo z 

" miriiQ " 


+ 




1 \ 


(niV\\Ti) + (nzV\\Tz 

miTiio L\ " / ' 


(13) 


In Eq. (13), {A) represents the ensemble average of a quantity A. In this study, we dehne 
an effective fluid velocity as (V||) = (V||j) + [{mznzo)/{Ywz)) which is the center of 
mass (GM) velocity in an equilibrium state of a combined ion-impurity system. and 

iZ^ 

nh are the contribution to the total parallel Reynolds stress from main ions and impurities. 


respectively. The last two terms in Eq. (13) represent the parallel turbulence acceleration 
whose characteristics have been studied elsewherd^lEs] 

In this study, we do not take the turbulence acceleration effect into account, focusing 
ourselves on the residual stress due to the symmetry breaking by the E x B shear. Also, we 
assume that the equilibrium flow velocities of main ions and impurities under consideration 
is equal, Vo||j = Vo\\z = Vo||- This assumption is strictly valid at the edge region of a tokamak 
plasma with sufficiently high collisionality. We keep this assumption in this study by taking 


notice of the validity of our analysis. Then, Eq. (13) is simplihed into the form, 

^(^l> 


dt 




(14) 
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where 11®® is the effective parallel Reynolds stress of the coupled ion-impurity system. One 
can write n®|f in a canonical form in the quasilinar theory, 



,,eff 


+ n 


eff 

res 


(15) 


with and n®f^ the effective momentum diffusivity and the residual stress dehned at the 
CM frame, respectively. Then, it is easy to derive expressions for y^® and n®f^ to obtain 


xf 


n®® = 

res 


i M z 

(16) 

n.e.+ 

(17) 


l + M 

where M = mznzo/miriio, (y|) is the turbulent ion (impurity) momentum diffusivity, 
and (nfes) is the residual stress whose radial gradient brings about the intrinsic torque. 
The physical origin for the appearance of M is obvious. Since we assume an equal flow 
velocity for main ions and impurities, the turbulence driven intrinsic torque is distributed 


between them according to their inertia. This effect is represented by M in Eqs. (16) and 


(1171). 

To complete the derivation, we need to evaluate the momentum diffusivities and parallel 
residual stress terms. It is straightforward to calculate them to obtain 


n*|| = {5veMu) = 


CsPs 




-Cn',, + ( h 

Kq, 


1 

1 + I 12 - 

car 


I 0 A 


^7 dr 


+ n* 


(18) 


nc = {6 veM\\z) = 


^sps 




k 






2kyUJrLezKz 

z\uj\^ 


\(t)k 


= -X 4 , 


z^^owz „z 

“T 


dr 


(19) 




n* 


n 


z 

res 


xl 


(^sPs 


CsPs 


CsPs 

Lue 

2 


V u |2 


ky (fc||) 


1 + 


2‘kyUjy^L^iIXi 


\UJ\ 


\<t^k\\ 




7- 


Z ky ^^11) f ^kyOJ^LgzKz 


P |a;| 


1 + 


Z|a;| 


M- 
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( 20 ) 

( 21 ) 

( 22 ) 


where 



























FIG. 5: The ratio of the residual stress (hires) to the diffusive momentum flux (If/)) as a 
function of fz- C®"'' (red) and (blue) are considered. Dotted and solid lines represent 
the cases when Lez = — 1 and 1, respectively. The quasi-linear momentum flux is 
calculated by integrating over all ky values giving rise to instability. Other plasma 
parameters being used are same as those of Fig. 2. 


Equations (15) to (22) are main results of the present work. They describe 11®® driven by 
ITG turbulence modihed by the presence of a single impurity species. 


The most important impurity effects on n®|[ are implicitly contained in Eqs. (20), (21), 


and (22). The characteristics of ITG turbulence, such as the growth rate, the mode shift 
giving rise to (fc||), and the fluctuation amplitude, change as impurity characteristics vary, as 


shown in Sec. III. Therefore, a complete study of impurity effects on 11®®^ requires calculations 
of actual value of |0 a:P by the variation of impurity characteristics. This is impossible in 
the quasi-linear theory. To overcome this limitation of the quasi-linear theory, we evaluate 
the ratio of n®f^ to the diffusive flux [11^ = —X^®(5Vo||/clr)] for a hxed equilibrium flow 
gradient. 


Teff 


n® 

D res 

ih — ~~ 

TTeff 


n 


eff 


(23) 


xfSVoi/dr' 

Then, R can be used as a quantitative measure, in the context of the quasi-linear theory, of 
the relative importance of the residual stress when impurity characteristics vary. Increase 
of R implies the grow of the relative importance of n^es compared to n^, and vice versa. 
Note that R does not involve |0 a:P term, which eliminates the difficulty of the quasi-linear 
theory in predicting the level of turbulence amplitude. 
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FIG. 6: The ratio of the residual stress (n^es) to the diffusive momentum flux (II^)) as a 
function of Lf,z = /in.. C«+ (red) and He^"*" (blue) are considered with a hxed value of 

jz = 0.3. A cyan dotted line represents the liresl^D value without impurities. Other 
plasma parameters being used are same as those of Fig. 5. 


Figure 1^ shows i? as a function of fz for C®’*' (red) and He^^ (blue) impurities. Plasma 
parameters being used to produce Fig. are same as those of Fig. 2. In calculating the 
radial momentum flux in Figs. 5 and 6, we keep all the unstable modes in Eqs. (20), (21), 


and (22) by integrating over all ky giving rise to positive growth rate. R shows a remarkably 
different trend depending on the alignment of an impurity prohle to the main ion density 
prohle. It increases when two prohles have a different sign (dotted lines) while it decreases 
when they have the same sign (solid lines). For example, the reduction in R as fz changes 
from 0 to 0.15 is 9.6% and 16.6% for He^^ and C®’*' cases, respectively. For a same absolute 
value of LeZ; the decrease rate of R when Lez > 0 is larger than the increase rate when 
LeZ < 0. Figure highlights the impact of the prohle alignment on R for a hxed value of 
fz = 0.3. The cyan dotted line in Fig. [^represents the R value when fz = 0. It shows a 
considerable reduction from its reference value (dotted line) as Lne/Lnz increases. 

To understand the results in Figs. |^and[^ we evaluate R for a hxed equilibrium velocity 
gradient and a ky value corresponding to the most unstable mode. Using the quasi-neutrality 
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condition, Eq. 0 becomes 


= 




l + M \L„ 

X ■! 1 + 

o; r 


E 


kyk\\ 


1 - 


fz 


Ki-fz[ K, 


LezKz 


\(t)k 


(24) 


where the relation M = mznzo/miriio = {p/Z) [/z/1 — fz] is used. Next, noting that 
= {x% + ^xf)/(l + M) = x^(j) from Eq. (21), one can easily show that 


{R)^ 


fell 


l + M \ ky/ I- fz 


1 I ‘^kyLJr J 
“T I lo 


\0J\ 


K,-fz{ K, 


LezKz 


(25) 


where {R)^ denotes an ensemble averaged value of R for a hxed ky. 


From the examination of Eq. (25), one can see that impurities affect R through four ways: 


(1) the mode shift which is contained in (k\\/ky'^, 

(2) impurity contents (fz), 

(3) stability change by the change of the sign of L^z, 

(4) the radial inhomogeneity of a impurity prohle represented by Kz, he., the direct 
impurity effect to the residual stress. 


All the hrst three effects turn out to be important in the determination of R. Among them, 
the mode shift, which is presented in Fig. 4, is found to be the most crucial factor determining 


the amount of residual stress by the variation of Lez- The last term in Eq. (25) represents a 
driving term to the residual stress solely from the free energy of impurities, he., the radial 


inhomogeneity of an impurity prohle. Positive L^z makes the last term of Eq. (25) negative 
since Lei > 0 in nominal tokamak discharges and < 0. Thus, this term can further reduce 
{R)f.^ combined with a strong reduction of |Xo| when Lez > 0. However, the effect of this 
term is marginal for cold impurity ions or high Z impurities like expected in machines 

with full tungsten PFC^. However, this term may play a role in {R) when fusion-born 
helium ashes are present because of their small charge number and natural tendency to 
align with main ions. Thus, results in this section might raise a potential concern to ITER 
operation where a signihcant amount of plasma rotation is envisioned to be produced via 
intrinsic rotation. 
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V. SUMMARY AND CONCLUSIONS 


In conclusion, this paper has shown that impurities may change the amount of intrinsic 
rotation by affecting the residual stress driven by turbulence. To make an analytic study, 
we employed a two fluid model for ions and impurities for ITG turbulence in sheared slab 
geometry. The ion and impurity velocities are assumed to be equal, which is strictly valid in 
the collisional regime. Then, we evaluated the momentum diffusivity and the residual stress 
in a coupled main ion-impurity system. The principal hndings of this paper are summarized 
as follows: 

• Characteristics of ITG modes, such as the growth rate or the mode structure, changes 
signihcantly by the variation of impurities contents and/or the alignment of an impu¬ 
rity prohle with that of main ions. Quasi-linear calculations show that these changes 
result in a considerable change of turbulence driven residual stress, hence the amount 
of intrinsic rotation. 

• The ratio of residual stress (n,.e<j) to the diffusive momentum flux (II^)), R = Ifres/II^), 
can be used as a quantitative measure of relative importance of If^es over 

• For a core-peaked impurity prohle shown in Fig. 1(a) schematically, R decreases due 
to the stabilization and a large increase of the mode shift of the unstable mode. For 
a wall-peaked impurity prohle shown in Fig. 1(b) schematically, R increases a little 
implying the predominance of If^es over II/). The degree of the enhancement of R 
when Lez < 0 is smaller than that of reduction when L^z < 0 for a given impurity 
content. 

As a possible implication of the present study, we pointed out that the amount of intrinsic 
rotation in burning plasma experiments might be reduced. This is due to the presence of 
He^^ ashes which is likely to be aligned with the main ion prohle. For example, approxi¬ 
mately ~ 10% of intrinsic rotation loss is expected if we assume the presence of He^’*' ions 
with fz = 0.15 well-aligned with a main ion prohle. Such a well-aligned impurity prohle is 
also possible in present day devices, as shown in recent experiments in KSTARP, due to the 
prevalence of turbulent impurity transport. We note that this may cause the decoupling of 
the ion temperature and the toroidal rotation prohles, which has been observed in recent 
KSTAR H-mode experiments^. When Lez > 0, the reduction of residual stress will be larger 
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than that of the diffusive momentum flux. Since the intrinsic torque is believed to be maxi¬ 
mized near the pedestal top region as shown in recent flux-driven gyrofluid simulationP^ES]^ 
the reduction of the residual stress due to impurities will result in the reduction of a net in¬ 
trinsic torque. This will reduce the total amount of toroidal rotation when an edge pedestal 
build up. To conhrm this scenario, measurements of a rotation prohle for a varying impurity 
contents is necessary. This is left as a future investigation. 

As an extension of the present study, we plan to perform a gyrokinetic analysis including 
impurities. Nonlinear gyrokinetic or gyrofluid simulations are also planned to elucidate the 
nonlinear physical process of intrinsic rotation generation and saturation in the presence 
of impurities. These will be outstanding research subjects which will be published in the 
future. 

ACKNOWLEDGMENTS 

This research was supported by the R&D Program through National Fusion Research Insti¬ 
tute (NFRI) funded by the Ministry of Science, ICT and Future Planning of the Republic 
of Korea (NFRI-EN1541-1), and by the same Ministry under the ITER technology R&D 
programme (INI504-1). 


17 




REFERENCES 


^D. G. Whyte, M. R.Wade, D. F. Finkenthal, K. H. Burrell, P. Monier-Garbet, B. W. Rice, 
D. P. Schissel, W. P. West, and R. D. Wood, Nucl. Fusion 38, 387 (1997). 

H. Lee (private communication, 2014). 

^A. Kallenbach, R. Dux, M. Mayer, R. Neu, T. Piitterich, V. Bobkov, J. G. Fuchs, T. 
Eich, L. Giannone, O. Gruber, A. Herrmann, L. D. Horton, G. F. Maggi, H. Meister, H. 
W. Muller, V. Rohde, A. Sips, A. Stabler, J. Stober, and ASDEX Upgrade Team, Nucl. 
Fusion 49, 045007 (2009). 

^F. Romanelli and JET EFDA Gontributors, Nucl. Fusion 53, 104002 (2013). 

®G. Federici, P. Andrew, P. Barabaschi, J. Brooks, R. Doerner, A. Geier, A. Herrmann, G. 
Janeschitz, K. Krieger, A. Kukushkin, A. Loarte, R. Neu, G. Saibene, M. Shimada, G. 
Strohmayer, and M. Sugihara, J. Nucl. Mater. 313-316, 11 (2003). 

®D. Reiter, H. Kever, G. H. Wolf, M. Baelmans, R. Behrisch and R. Schneider, Plasma 
Phys. Gontrol. Fusion, 33, 1579 (1991). 

^A. R. Polevoi, M. Shimada, M. Sugihara, Yu. L. Igitkhanov, V. S. Mukhovatov, A. S. 
Kukushkin, S. Yu. Medvedev, A. V. Zvonkov, and A.A. Ivanov Nucl. Fusion 45, 1451 
(2005). 

®S. P. Hirshman and D. J. Sigmar, Nucl. Fusion 21, 1079 (1981). 

W. B. Kim, P. H. Diamond, and, R. J. Groebner, Phys. Fluids B 3, 2050 (1991). 

^°R. Dux and A. G. Peeters, Nucl. Fusion 40, 1721 (2000). 

Garbet, N. Dubuit, E. Asp, Y. Sarazin, G. Bourdelle, P. Ghendrih, and G. T. Hoang, 
Phys. Plasmas 12, 082511 (2005). 

^^G. Angioni and A. G. Peeters, Phys. Rev. Lett. 96, 095003 (2006). 

^^B. Goppi, H. P. Furth, M. N. Rosenbluth, and R. Z. Sagdeev, Phys. Rev. Lett. 17, 377 
(1966). 

Q. Dong, W. Horton, and W. Dorland, Phys. Plasmas 1, 3635 (1994). 

Q. Dong, and W. Horton, Phys. Plasmas 2, 3412 (1995). 

^®X. Y. Fu, J. Q. Dong, W. Horton, G. T. Ying, and G. J. Liu, Phys. Plasmas 4, 588 (1997). 
^^T. Fiilop and J. Weiland, Phys. Plasmas 13, 112504 (2006). 

^®G. Angioni, L. Garraro, T. Dannert, N. Dubuit, R. Dux, G. Fuchs, X. Garbet, L. Garzotti, 
G. Giroud, R. Guirlet, F. Jenko, O. J. W. F. Kardaun, L. Lauro-Taroni, P. Mantica, M. 


18 



Maslov, V. Naulin, R. Neu, A. G. Peeters, G. Pereverzev, M. E. Puiatti, T. Piitterich, 
J. Stober, M. Valovic, M. Valisa, H. Weisen, A. Zabolotsky, ASDEX Upgrade Team, and 
JET EEDA Gontributors, Phys. Plasmas 14, 055905 (2007). 

Nordman, R. Singh, T. Eiilop, L.-G. Eriksson, R. Dumont, J. Anderson, P. Kaw, P. 
Strand, M. Tokar, and J. Weiland, Phys. Plasmas 15, 042316 (2008). 

^°T. Eiilop and H. Nordman, Phys. Plasmas 16, 032306 (2009. 

M. Tang, R. B. White, and P. N. Guzdar, Phys. Fluids 23, 167 (1980). 

^^P. H. Diamond, G. J. McDevitt, O . D. Giircan, T. S. Hahm, W. X. Wang, E. S. Yoon, I. 
Holod, Z. Lin, V. Naulin, and R. Singh, Nucl. Fusion 49 , 045002 (2009). 

^^P. H. Diamond, G. J. McDevitt, 6 . D. Giircan, T. S. Hahm, and V. Naulin, Phys. Plasmas 
15 , 012303 (2008). 

^^6. D. Giircan, P. H. Diamond, T. S. Hahm, and R. Singh, Phys. Plasmas 14 , 042306 
(2007). 

D. Giircan, P. H. Diamond, P. Hennequin, G. J. McDevitt, X. Garbet, and G. Bourdelle, 
Phys. Plasmas 17 , 112309 (2010). 

^®G. J. McDevitt, P. H. Diamond, 6. D. Giircan, and T. S. Hahm, Phys. Rev. Lett, 103 , 
205003 (2009). 

^^Rameswar Singh, Rajaraman Ganesh, Raghvendra Singh, Predhiman Kaw, and Abhijit 
Sen, Nucl. Fusion 51 , 013002 (2011). 

^®Rameswar Singh, R. Singh, Hogun Jhang and P. H. Diamond, Phys. Plasmas 21 , 012302 
(2014). 

^®R. Balescu, Aspects of anomalous transport in plasmas, (TOP Publishing Ltd., Bristol, UK, 
2005), p. 55. 

^°P. Gohil, K. H. Burrell, and T. N. Garlstrom, Nucl. Fusion 38 , 93 (1998). 

^^N. G. Hawkes, D. V. Bartlett, D. J. Gampbell, N. Deliyanakis, R. M. Gianella, P. J. Lomas, 
N. J. Peakcock, L. Porte, A. Rookes, and P. R. Thomas, Plasma Phys. Gontrol. Fusion 
38 , 1261 (1996). 

^^L. Wang and P.H. Diamond, Phys. Rev. Lett. 110 , 265006 (2013). 

^^X. Garbet, D. Esteve, Y. Sarazin, J. Abiteboul, G. Bourdelle, G. Dif-Pradalier, P. Ghen- 
drih, V. Grandgirard, G. Latu, and A. Smolyakov, Phys. Plasmas 20 , 072502 (2013). 

34Won-Ha Ko, H. Lee, P.H. Diamond, S.H. Ko, J.M. Kwon, L. Terzolo, S. Oh, K. Ida, Y.M. 
Jeon, K.D.Lee, J.H. Lee, S.W. Yoon, Y.S. Bae, W.G. Kim, Y.K. Oh, and J.G. Kwak, in 


19 



Proceedings of the 24 th IAEA Fusion Energy Conference, San Diego, USA, October 8-13, 
2012 (International Atomic Energy Agency, Vienna, 2012), PD/P8-22. 

^^Hogun Jhang, S. S. Kim, and P. H. Diamond, Jour. Kor. Phys. Soc. 61, 55 (2012). 

^®S. Tokunaga, Hogun Jhang, S. S. Kim, and P. H. Diamond, Phys. Plasmas 19, 092303 

( 2012 ). 


20 



